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Using X-ray photon correlation spectroscopy, we studied the dynamics in the nematic phase of a nanorod 
suspension. The collective diffusion coefficient in the plane perpendicular to the director varies sharply with 
the wave vector. Combining the structure factor and the diffusion coefficient we find that the hydrodynamic 
function of the phase decreases by more than a factor of ten when going from length scales comparable to the 
inter-particle distance towards larger values. Thus, the collective dynamics of the nematic phase experiences 
strong and scale-dependent slowing down, in contrast with isotropic suspensions of slender rods or of spherical 
particles. 
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Nematics are the simplest example of a phase with no 
positional order, but still exhibiting orientational order 
(and hence anisotropy). This combination endows them 
with remarkable qualities: although fluid, they have elas¬ 
tic properties and, consequently, long-lived fluctuations. 
A great deal is known about the large-scale behaviour of 
nematic systems, which is well-described by a generalized 
hydrodynamic modeli. This theoretical description was 
confirmed (and refined) using a wealth of experimental 
techniques. The method of choice for studying nema- 
todynamics is dynamic light scattering (DLS), which is 
sensitive to the relaxation of nematic fluctuations on mi¬ 
cron scales^i^. 

On the other hand, there is much less data on the 
short-range dynamics of the nematic phase, convering 
length scales comparable to the inter-particle distance. 
In this limit, the continuous medium model is bound 
to break down, and more microscopic considerations 
must be taken into account. Since this is the scale at 
which interaction between particles defines the structure 
of the system, understanding the dynamics is essential 
for building a complete picture of the phase. A consid¬ 
erable body of theoretical and numerical work exists^i^, 
but there is very little experimental data, mainly due to 
the lack of adapted techniques (due to the typical par¬ 
ticle size, this range of scattering vectors is out of reach 
for DLS.) Alternative methods can be used, such as in¬ 
elastic neutron scattering, which is however limited to 
sub-microsecond dynamics (too fast for cooperative pro¬ 
cesses). Spin relaxation has also been employed, but it 
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lacks the required space resolution and the conclusions 
are indirect. 

X-ray scattering techniques are suitable for exploring 
these distances, but until recently they were only able 
to draw a static picture of the system. This situation is 
changing due to the progress of X-ray photon correlation 
spectroscopy (XPCS), opening up the time dimension. 
However, due to inherent technical difficulties, the ex¬ 
perimental systems must fulfill very stringent conditions, 
such as slow relaxation rates and high scattering con¬ 
trast. Some experiments have already been performed 
on nematics using XPCS, but they were only concerned 
with very slow relaxation in a gel phased or with capillary 
surface waves^. 

Essential information on the physics of multi-particle 
systems is contained in the structure factor S{q) and the 
collective diffusion coefficient D{q)^. Based on very gen¬ 
eral thermodynamic arguments, these parameters are re¬ 
lated by D{q) ^ l/*S'(g) (“de Gennes narrowing”—). In 
colloidal suspensions, a more refined treatment must take 
into account the hydrodynamic interactions, which fur¬ 
ther modulate this dependence. 

In this letter, we study a fluid nematic phase of 
goethite (a-FeOOH) nanorods (with moderate aspect ra¬ 
tio) and determine the hydrodynamic function over a q- 
range corresponding to length scales comparable to the 
inter-particle distances. Unexpectedly, we find that for 
wave vectors q < ^max (below the maximum of the struc¬ 
ture factor), the dynamics of the system slows down con¬ 
siderably. This result is in stark contrast with isotropic 
suspensions of slender rods, where no hydrodynamic ef¬ 
fect is observedi^. Furthermore, the effect is much 
stronger than the variation of the hydrodynamic func¬ 
tion in suspensions of colloidal spheres^i, emphasizing 
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the role of particle anisotropy and of the nematic order. 
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FIG. 1. a) SAXS pattern of an aligned nematic sample with 
(pg — 6.7% at 228 K. The arrow indicates the direction of the 
nematic director n. b) Structure factor along the q± direction 
calculated from the SAXS data in a). 

Goethite nanoparticles were synthesized according to 
a well-established procedure and dispersed in wate r^^’^^ . 
The particles are lath-shaped; transmission electron mi¬ 
croscopy images show that they are 4000 ± 1000 A long 
and 330 ± 110 A wide ( mean and standard deviation 
obtained from a log-normal fit). Moreover, the particle 
cross-section is usually anisotropic by a factor of 2.5^^. 
In the following, we therefore consider that the particles 
rotate freely around their main axis, as expected in a uni¬ 
axial nematic phase. In order to slow down their dynam¬ 
ics so that it may be conveniently studied by XPCS, the 
particles were resuspended in an 80 wt.% propane-1,3- 
diol in water mixture. Although this mixture has lower 
dielectric permittivity than water, the particles still bear 
some positive electrical surface charge density, which en¬ 
sures the colloidal stability of the suspension. In the 
first approximation, the effect of electrostatic interactions 
can be described by introducing an effective radius larger 
than the bare one, which tends to reduce the apparent 
particle aspect ratio^^. For example, previous studies of 
another goethite suspension in pure water have shown 
that the ratio of effective diameter to bare diameter is 
about 1.6i^. With the solvent mixture used here, this 
electrostatic correction factor is expected to be smaller. 


Therefore, in the following, the particles will be consid¬ 
ered as cylindrical rods, 4000 A long and with an aspect 
ratio of about 10. (The interpretation presented below 
does not critically depend on the precise value of the ef¬ 
fective aspect ratio.) At volume fractions (j)g < 4.2% the 
suspensions are isotropic, whereas at (fig > 6.7% they 
form a nematic liquid-crystalline phase that aligns in low 
magnetic field, with its nematic director parallel to the 
fieldi^. Samples of different volume fractions were held 
in optical flat glass capillaries, 50 /j,m thick (VitroCom, 
NJ, USA) and placed in a vacuum chamber. The tem¬ 
perature was lowered to 228 K, where the viscosity of the 
solvent is of the order of 1000 mPa.s. The nematic phase 
was aligned with a 150 mT field that was then removed 
for the actual measurements. 

The small-angle X-ray scattering (SAXS) and XPCS 
measurements were performed at the TROIKA beam line 
IDIOA of the ESRF with an X-ray energy of 8 keV (A = 
1.55 A) selected by a single-bounce Si(lll) monochro¬ 
mator, in the uniform filling mode of the storage ring. A 
(partially) coherent beam is obtained by inserting a 10 
/im pinhole aperture a few centimetres upstream of the 
sample. 

For the XPCS measurements, we used a two- 
dimensional (2D) Maxipix detector consisting of 256 x 256 
square pixels (55 /im in size), and the intensity autocor¬ 
relation functions were calculated by ensemble averaging 
over equivalent pixels^i. In the nematic phase, the pixels 
averaged were restricted to a narrow slice perpendicu¬ 
lar to the nematic director for qj_ (see Figure [T] for an 
illustration and Figure [2] for the results), or parallel to 
the nematic director for qy. In the isotropic phase, all 
pixels at the same scattering vector modulus were aver¬ 
aged. Some measurements were also performed using a 
point detector (an avalanche photodiode) connected to 
an external digital correlator (FlexOlD-08). 

A typical SAXS pattern of an aligned nematic sam¬ 
ple is shown in Figure [1^. The static structure factor 
S{q±) in the q± direction (Figure [T}d) was obtained by 
dividing the scattered intensity by the form factor mea¬ 
sured independently on a dilute solution. S{q±) dis¬ 
plays a well-defined interaction peak, at a value (/max = 
8.6 X 10“^ A due to the liquid-like positional short- 
range order of the nanorods in the plane perpendicular 
to the director. 

The various dispersion relations are shown in Figure 
O The relaxation rate in the nematic phase, perpendic¬ 
ular to the director {N , T, solid triangles) is linear in 
roughly up to the position of the structure peak. At 
higher q±, the slope increases abruptly, before approach¬ 
ing a final linear regime. In contrast, the behaviour along 
the director (N , ||, open diamonds) is linear over the ac¬ 
cessible range, which is limited by the rapid fall-off of 
the intensity in this direction (see Figure [TJi). In the 
isotropic phase (with (jig = 2%), the dispersion relation 
is also linear over the entire range. 

Let us define the collective diffusion coefficient 
DN,±{q±) = ^{q±)/Q±, shown in Figure |3 Jd). For ref- 
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FIG. 2. Autocorrelation functions for a nematic sample with 
= 6.7% at different scattering vectors q±. The solid line is 
the ht with a stretched exponential (the stretching exponent 
is about 0.6 for all the curves). 



FIG. 3. Dispersion relations for a nematic {(pg = 6.7%) sam¬ 
ple, with q parallel to the nematic director (o), and with q 
perpendicular to the nematic director (a). The dispersion 
relation of an isotropic, {cfig = 2%) sample is shown for com¬ 
parison (x). 


there is no physical reason for adopting it. 



erence, panel a) shows the structure factor (Figure [TJr). 

Two regimes (above and below the interaction peak of 
the structure factor) can be clearly distinguished. Be¬ 
tween them, the diffusion coefficient jumps by more than 
a factor of 3. To quantify this variation, we fit DN ±{q±) 
with a sigmoidal function (solid line in Figure IUd): 


D{q) — Umin + 


i + »p(-^) 


( 1 ) 


with parameters: Umin = 1.1 x 10 and Umax = 
3.4 X 10-^®mVs while qi /2 = 9.6 x and Aq = 
0.45 X 10“^ A“^. The greyed areas around the fit are the 
±(T prediction bands, quantifying the data scatter (about 
68 % of the experimental points should fall within this 
area). The fit function dl]) is only chosen for convenience; 


FIG. 4. a) Static structure factor S{q±). b) Diffusion coef¬ 
ficient DN,±{q±)- The symbols correspond to two different 
measurements performed on the same sample. The solid line 
and the greyed area are the fit by a sigmoidal function and 
the ±CT prediction bands (see text), c) Hydrodynamic func¬ 
tion obtained using Equation [21 The data are plotted against 
a common q± axis. 

The decrease of DN^±{q±) at low wave vectors is even 
more striking if we recall that, for q± vectors below the 
peak, where the structure factor decreases, the diffusion 
coefficient should increase\ We can restate this result 
more precisely in terms of the hydrodynamic function, 
H{q±), defined hy^: 

^ £N^Ms(q^) ( 2 ) 

and shown in Figured!:;). 
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Care must be exercised when determining the constant 
factor -Do,_Lj which is the diffusion coefficient of the par¬ 
ticles at infinite dilution. The _L subscript specifies that 
one should consider the diffusion in the direction perpen¬ 
dicular to the main axis of the particle (in the nematic 
phase this axis is parallel to the director, and we are 
now concerned with diffusion perpendicular to the direc¬ 
tor). We start by measuring the (orientationally aver¬ 
aged) diffusion coefficient in the isotropic phase, Dj. For 
4>g = 2%, the data is shown in Figure [3] (crosses), along 
with a linear fit yielding Di{2%) = 6.4 x 10“^®m^s“^. 
Similar measurements at (j)g = 0.5% (data not shown) 
give Dj{0.5%) = 8.0 x 10“^® m^ s“^. Correcting these re¬ 
sults for the intrinsic viscosity of nanorod suspensions^^ 
leads to an infinite dilution value Do = Di{0) = 8.5 x 
10“^®m^s“^. In isotropic solution, Dq is a geometric 
combination of the diffusion coefficients along and per¬ 
pendicular to the major axis: SDq = Do,|| + 2Do,± and, 
for nanorods with an aspect ratio of l(), Dq ,±/Do = 
0.88^^, resulting in Z3 o,-L = 7.5 x 10“^® m^ s“^. 

At high q±, H{q) reaches a value lower than one 
{H{oo) ~ 0.4). For spherical particles, this limit is de¬ 
scribed by2i: H{oo) = rjo/rj, where ry is the viscosity of 
the suspension and rjo that of the solvent. Transposing 
this formula to our (non-spherical) system yields fairly 
good agreement^^. The main feature of H {q) is however 
the decrease at lower wave vector, starting around (/max 
and clearly visible in Figure |4j:). At the lowest accessi¬ 
ble wave vector H{q) ~ 0.04, ten times smaller than the 
maximum value H{oo). 

Hydrodynamic slowing down of the collective relax¬ 
ation is also encountered in suspensions of spherical par¬ 
ticles, but its amplitude is much lower; for a volume 
fraction (j) = 9% (higher than in our nematic phase), 
H{0)/H{oo) ~ l/3ii. Experimental and theoretical re¬ 
sults for even higher volume fractions of spheres (both in 
the low- and high-salt concentration regimes) yield more 
modest decreases for the collective diffusiort^^. We con¬ 
clude that the behaviour of H{q) at low wave vectors in 
the nematic phase is very different from that in sphere 
suspensions. 

A possible explanation for this difference is that, if the 
major axis of the particles is much longer than the typi¬ 
cal distance between them, (in the limit of very large or¬ 
der parameter and aspect ratio)the nematic phase should 
behave like a two-dimensional (2D) system^!, where hy¬ 
drodynamic interactions are stronger than for a three- 
dimensional systent^^. This explanation is in qualitative 
agreement with recent simulations of 2D colloidal suspen¬ 
sions, where the hydrodynamic interactions slow down 
the collective diffusion coefficient^^. In particular, for the 
volume fraction used in our study these authors find that 
the hydrodynamic interactions slow down the collective 
diffusion by at least a factor of four (Ref. 26, Figure 2b). 
However, more theoretical or numerical results would be 
needed, in particular concerning the length-scale depen¬ 
dence of the diffusion coefficient D{q), to understand this 
pronounced slowing down. 
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